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A study of the finite-size effect is carried out for spectra of both ground-state 
and excited-state baryons in quenched lattice QCD using Wilson fermions. Our 
simulations are performed at (3 = Q/g 2 = 6.2 with three different spatial sizes, 
La ~ 1.6, 2.2 and 3.2 fm. It is found that the physical lattice size larger than 3 
fm is required for A states in all spin-parity ( J p = 1/2 ,3/2^) channels and also 
negative-parity nucleon ( N *) state ( J p = l/2 _ ) even in the rather heavy quark mass 
region (M n ~ 1.0 GeV). We also find a peculiar behavior of the finite-size effect on 
the hyperfine mass splittings between the nucleon and the A in both parity channels. 

PACS numbers: 11.15.Ha, 12.38.-t 12.38.Gc 


I. INTRODUCTION 


In the past several years, many lattice studies on excited baryon spectroscopy have been 
made within the quenched approximation after the pioneering attempts jl|. Although these 
calculations gave results of negative-parity baryon masses in gross agreement with exper¬ 
iments (2, 13, 1H E, 3 H ft, 10 • the higher precision study requires accurate controls of 
systematic errors, which may arise from the quenched approximation, chiral extrapolation, 
nonzero lattice spacing and finite-size effect. The first and second errors might appear most 
seriously in spectra of the excited states and should be entangled with each other due to 
possible decay processes. This complicated issue, however, ought to be postponed since 
the current dynamical simulations are still limited even for the ground-state baryons. It is 
worth emphasizing that the quenched lattice QCD is still useful to investigate several long 
standing puzzles in excited baryon spectroscopy, such as the level ordering problem between 
the negative-parity nucleon 1V*(1535) and the Roper resonance A^'(1440), the structure of 


the A(1405) and so on 11 


Instead of using small lattice spacing, many calculations have adopted sophisticated lat¬ 
tice discretization schemes for fermions, the improved fermion actions 3. 0[ E, Q or chiral 
fermions Jjj, 7, E, 3], hr excited baryon spectroscopy to reduce the cutoff effect and improve 
the chiral behavior. However, there is not so much attention paid to the finite-size effect. 
Indeed, most of all simulations are performed with spatial lattice sizes La ~ 1.6 — 2.2 fm. 
This is because the finite-size study has been hardly done due to large computational costs. 
In this study, we utilize the standard Wilson fermion action, the computational requirements 
of which are relatively cheap. Therefore, we systematically perform the finite-size study for 
excited baryon spectroscopy. 

For the finite-size effect, after the pioneering analytic consideration 42j|, where the finite- 
size behavior of hadron masses obeys the exponential form M ^ — Ml ~ L -1 exp(— L/Lq), 
it has been found that this behavior is rather described by a power law ~ L~ n with n = 
2-30. This power-law behavior can be interpreted by the “wave function” squeezed on 
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a finite lattice from the phenomenological point of view. The dependence of the f ini te-size 
effect on the spatial boundary condition which is found earlier [III,El, can also be explained 
by this interpretation [l 6 |. Ref. [u| reported that the lattice size La > 2.5 fm is required 
for the nucleon ground-state (N) in quenched lattice QCD. 

In the non-relativistic quark model with a harmonics-oscillator potential H = ^p 2 + 
^mu 2 r 2 , the energy level and the root mean square radius can be estimated as E N = 


ku(N + |) and r rms = y + §), respectively. The negative-parity nucleon (. N *) is 

identified as N = 1 state, which corresponds to the first orbital excitation. It turns out 
that a ratio of the root mean square radii between the N* and the N states is estimated as 
'i/'i = vT's 1.3. This means that the finite-size effect stemming from the hadron 
size squeezed is expected to become much serious for the excited-state baryons rather than 
the ground-state baryon. According to a crude estimate with this simple model, the lattice 
size La > 2.5 x 1.3 ~ 3.2 fm might be required for the excited baryon spectroscopy. 

The organization of our paper is as follows. In Sec 01 we first review the basic idea of the 
parity projection in lattice simulations. In recent years, the parity projection technique is 
greatly appreciated in determining the parity of the pentaquark state in lattice QCD [171, 18]. 
Unfortunately, however, there are some initial confusions about the parity assignment |17| . 
Thus, we describe details of the precise parit y p rojection, which is always adopted by one 
of authors in excited baryon spectroscopy [3. ful 19] - For the spectroscopy of the spin-3/2 
state, he. the A baryon, appropriate spin projections are also described. Sec. II I II gives 
details of our Monte Carlo simulations and results of the f ini te-size study for various baryon 
states, such as the nucleon states in both parity channel ( J p = 1 /2 ± ) and the A states in 
all spin-parity channels ( J p = l/2 ± ,3/2 ± ). Then, we discuss about mass splittings between 
hyperhne partners in both parity channels. Finally, in Sec. I1VI we summarize the present 
work. 


II. GENERAL ANALYTIC FRAMEWORK 
A. Parity projection 

In general, the local baryon operator should couple to both parity states. The local 
baryon operator is defined by the trilinear composite operator, which is composed of a local 
diquark operator and a spectator-like quark held 0 as 

Ob{x) = £abc (qli{x)CTq b j(x)) T'q C}k (x) (1) 

where T and T' stand for possible 16 Dirac matrices and C is the charge conjugation matrix. 
The superscript T denotes transpose. Indices abc and ijk have meanings as color and flavor 
respectively. The intrinsic parity of this operator is defined by the parity transformation of 
internal quark fields Vq(x,t)V^ = + 74 q(—x,t) as 

VO$\x,t)V' = r i 'Y 40 $ ) (-x,t), (2) 

where p denotes the intrinsic parity of Ojf . However, due to the simple relation between the 
positive- and negative-parity operators: 0^\x) = 75 Og \x) for the local baryon operator , 
the resulting two-point correlation functions are also related with each other [ 2 | as 

(OloSt+JoSt+JIO) = -75<0|0<r+)oS 3 “ ) ( !/ )|0)7 5 . 


(3) 
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This means that the two-point correlation function constructed from the local baryon oper¬ 
ator has overlap with both parity states. Next, we consider the spin-1/2 baryon case as a 
typical example for the precise parity projection. The zero-momentum two-point function 
is given by the sum over all spatial coordinates x, 


G { S\t) = X>|T{0<?Ri)og , (0,0)}|0} ■ 

X 


This two-point correlation has the asymptotic form 


21 


G^\t) = ^ (1 + sgn(t) 74 ) e~ M ^\ - (1 - sgn(f) 74 ) e~ M ~^ 


(4) 

(5) 


at large Euclidean time t. Here M v and M_ v denote masses of the lowest-lying state in 
each parity channel. The amplitude A v (A_, ; ) is defined as (0\O B ] \B, rj) = y/A^u^\ Here, 
it is worth noting that anti-particle contribution of the opposite parity state is propagating 
forward in time. Thus, the +/— parity eigenstates in the forward propagating contribution 
can be obtained by choosing the appropriate projection = (1 d= 77 74 )/ 2 , which is given in 
reference to the intrinsic parity of the utilized operator, rj. This procedure is accomplished 
by taking the trace of the correlator over spinor with the relevant projection operator as 
t-Tr {P±^G B \t)} where the factor of 1/4 is our choice of normalization. 

The lattice simulation is, however, performed on a lattice with finite extent T in the time 
direction and (anti-)periodic boundary conditions. The definition of the forward propagating 
contribution is slightly ambiguous since the backward propagating can wrap the whole time 
range around the time boundary. Hereafter, for a simplicity, we take 77 = +. Under the 
(anti-)periodic conditions, the resulting correlation functions may be given by 


OO 

Gg b ' c '(t) = E G B {t + nT), (6) 

n=—oo 

00 

G‘ B p bC W = E (-)"G B (t + nT), (7) 

n =—00 

which certainly satisfy G p B h ' c '(t ) = G^ h c ' (t + T) and G a B p ' h ' c '(t) = —G^ h ' c '(t + T) respectively. 
For the region 0 <t<T, those correlation functions are expressed in the following forms: 


G p B hc (t) = E^+e 


—nM+T 


n =0 


-J2 A - e 


-nM-T 


n =0 


cy-h-'-w = E^ 

n=0 


j—nM+T 


-J2 A - e 


-nM-T 



2 

T 

'1 - 74 


2 

'1 +74„_ 


2 

T 

1-74 


+ 


e-M.t + 


-M+(T-t) 




-M+(T-t) 


-( T-t ) 


( 8 ) 


(9) 


n=0 L - 1 

For details of above calculation, see Appendix El In the large T limit (M + T 1 and 
M-T 1), all n/0 terms are negligible in a summation with respect to n, because of the 
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factor e nM ± T . Thus, it turns out to be the form 


G 


p.b.c./a.p.b.c. 
B 


it) « A, 


1+74 M+ t ± l _ 74 -M+(T-t) 

2 2 


-A. 


\ _74 -Af_t _j_ 1 + 74 -M_(T-t) 


( 10 ) 


which has only a contribution from t/ie first wrap-round effect. Clearly, even after the 
“parity projection’’ as lTr{P±G^ b ' c ^ a ' p ' b ' c '(t)}, there remain unwanted contaminations from 
the backward propagating contributions due to the first wrap-round effect. However, one 
can easily find that the linear combination of two correlation functions with periodic and 
anti-periodic boundary conditions may prevent the first wrap-round effect: 

a B (t) = i{G£ b -°'(«) + G %^'(«)} » (i±+) A + e- M -< _ (Wli) A.e~ M -‘. (11) 


Therefore, the parity projection are really accomplished through the operation 
iTr {P ± G B {t)}. 

To achieve this in an efficient manner, we adopt a linear combination of the quark propaga¬ 
tors with periodic and anti-periodic boundary conditions to construct the baryon two-point 
correlation function [i]. The linear combination in the quark level automatically realizes 
the linear combination in the hadronic level, as shown below. First, the baryon two-point 
correlation function G B {t) can be schematically written by 


Gs(t) = £ tr c {S'(a;, 0) • S(x, 0) • S(x, 0)} (12) 

X 

with the quark propagator S(x, 0) where x = (x,t) and 0 = (0,0). tr c denotes a trace over 
color indices. We define the averaged quark propagator 

S(x, 0) = i {S P {x, 0) + S'ap^O)) , (13) 

where Sp(x, 0) and Sap(x, 0) are subject to the periodic and anti-periodic boundary con¬ 
ditions in the time direction such as Sp(x,t] 0) = Sp(x,t + T; 0) and S'ap(x, t; 0) = 
—SAp(x,t + T; 0). By inserting this linear combination into Ea. (1121) . the baryon two-point 
function is written as 

G B (t) = tr c {(Sp(x, 0) + S A p(x, 0)) • (S P (x, 0) + Sap(z, 0)) ■ (S P (x, 0) + S A p(x, 0))} 

X 

= l{Gg b -'(*)+Gy bc (i)}, (14) 


where 

G p B h c ' (t) = i ^tr c {5 P (a;, 0) • S P (x, 0) • S P (x, 0) + S P (x, 0) • S A p(x, 0 ) • -S'ap^O) 

X 

+S A p(x, 0) • Sp(x, 0) • S A p(x, 0) + Sap(x, 0) • S A p(x, 0) • S P (x, 0)} (15) 
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and 


G a B p hc '(t) = ^ tr c {S AP (x, 0) • S P {x, 0) • S P (x, 0) + S P (x, 0) • S AP (x, 0) • S P (x, 0) 

X 

+S P (x, 0) • S P (x, 0) • S A p (ic, 0) + S AP (x, 0) • S AP (x, 0) • S AP (x, 0)} . (16) 

Gg b ' c '(i) has even numbers of S AP (x, 0) and then are satisfied with G p B h ' c '(t) = G^ h ' c '(t + T). 
On the other hand, G^ pbc '(t) includes odd numbers of S AP (x, 0) and is totally subject to the 
anti-periodic boundary condition, Gg P ' hc '(t) = —G^ p b c '(f + T ). We confirm that a linear 
combination of the quark propagators with periodic and anti-periodic boundary conditions 
makes a cancellation for the primal reflection from the time boundary in the baryon two- 
point function. This is also true for the meson two-point function as same as the baryon 
case. 


B. Spin projection 

In this subsection, we describe the spin projection, which is essential to deal with spin-3/2 
baryon, namely the A state. We choose a interpolating operator for the A, more specifically 
A ++ , as a simple trilinear composite operator: 

0^{x) = e abc (v%(x)C'y IJ ,u b (x)) u c (x), (17) 

where u(x) represents the up quark field. This operator has the structure of the Rarita- 
Schwinger spinor (“vector-spinor”) so that the two-point function constructed from this 
operator can couple to both spin-3/2 and spin-1/2 states 22], We consider the two-point 
function at zero spatial momentum, which is given by 

= X>|T{0£(f,«)0?(0,0)}|0). (18) 

X 

Note that above A correlator can couple to the spin-3/2 state only if neither of /j,,u are 
temporal indices. For the spatial Lorentz indices, i , j = 1,2,3, Eq. m is expressed by the 
orthogonal sum of spin-3/2 and spin-1/2 components: 

G?j(t) = (% - G$ /2 (t) + pOiGt/M, (19) 

with appropriate spin projection operators [2^]. Respective spin parts, G^ 2 (t) and G^ 2 (t), 
possess both positive- and negative-parity contributions as same as Eq. ©• Some of excited 
A states such as the spin-1/2 A states can be accessed by the correlator Gy 2 (t). These spin 
projected correlators are given by 

Gf n {t) = 5 G g(t) - W 7 atGiM, (20) 

k 

Gf /2 (i) = J>7*G£(t), 

k 


( 21 ) 
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where an index k should be summed over all spatial directions, but any specific choice is 
available for an index i. In this paper, we calculate all three direction for the spatial index 
i and take an average of them to get the possible reduction of statistical errors. 

According to the previous subsection, we perform the appropriate parity projection to 
apply to Gf, 2 (t) and Gy 2 (t) as well. Then, we finally access to the four different spin-parity 
states of the A baryon, which correspond to the spin-3/2 positive-parity state (A 3 / 2 ), the 
spin-3/2 negative-parity state (Ag/ 2 ), the spin-1/2 positive-parity state (Ax/ 2 ) and the spin- 
1/2 negative-parity state (A*y 2 ). Of course, the A 3 / 2 is the ground state of the A baryon. 
All others are the excited A baryons. 


III. NUMERICAL RESULTS 

We generate ensembles of the quenched QCD configurations with the standard single- 
plaquette action at /3 = 6 /g 2 = 6.2 with three different lattice sizes, L 3 x T = 24 3 x 48, 
32 3 x 48 and 48 3 x 48, and compute the quark propagators by using Wilson fermion action 
at several values of the hopping parameter k, which cover the range M n ss 0.6 — 1.2 GeV. 
Details of simulations are summarized in Tables |T| and HD 

For the update algorithm, we utilize the Metropolis algorithm with 20 hits at each link 
update. The first 10000 sweeps are discarded for thermalization. The 0(200-300) gauge 
ensembles in each simulation are separated by 1200 (L=48), 800 (L— 32) and 600 (L=24) 
sweeps. For the matrix inversion, we use BiCGStab algorithm [ 24 J and adopt the convergence 
condition |r| < 10~ 8 for the residues. We calculate a simple point-point quark propagator 
with a source location at f src = 6. To perform the precise parity projection, we adopt a 
procedure to take an average of two quark propagators which are subject to periodic and 
anti-periodic boundary conditions in time, as described in Sec. [H] We use the conventional 
interpolating operators, £ abciu^C^db)u c for the nucleon and £ a bc(u^C'y fJ ,Ui ) )u c for the A 
respectively. 

All calculations were done on a Hitachi SR8000 parallel computer at KEK - High Energy 
Accelerator Research Organization, using the extended code based on the Lattice Tool Kit 
(LTK) |25]. I 11 this paper, we analyze the baryon-mass spectra in either parity channels 
with the conventional single exponential fit. 


A. Spatial lattice-size dependence 


We first discuss the finite-size effect on masses of all measured baryons. It is important 
to investigate how large size of physical lattice is required to neglect the finite-size effect on 
exited baryon spectroscopy. For this purpose, we perform numerical simulations at three 
different lattice sizes, L 3 x T = 24 3 x 48, 32 3 x 48 and 48 3 x 48 for k = 0.1506 and 0.1497. 
Quenched f3 = 6.2 corresponds to a lattice cutoff of cU 1 = 2.913 GeV, which is set by the 
Sommer scale 0. Thus, the spatial extents in our study correspond to La ~ 1.6, 2.2 and 
3.2 fm in the physical unit. We will calculate the fini te-size corrections to the baryon mass 
on each finite volume by comparison with values evaluated in the infinite volume limit. 

We utilize the phenomenological powe- law formula to take the infinite volume limit for 
the observed masses as 


uMl = aihoo + cL 


(22) 
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This power-law behavior can be interpreted by the hadron “wave function” squeezed on a 
finite lattice from the phenomenological point of view 0 As for the quenched simulation, 
Kef. Ui reported that the range of n — 1 — 2 is preferable for the ground state of the nucleon. 
In this paper, we use the three different powers, n — 1,2 and 3, and then determine which 
power-law behavior is most favorable for fitting data. 

In Figs ITHHl the lattice-size dependence of each baryon mass (N, N *, A 3 / 2 , AW, Ax/ 2 , 
AW) is shown for two hopping parameters (k = 0.1497 and 0.1506), which correspond to 
the relatively heavier quark masses. All data included in Figs. 171131 are summarized in Table 
inn The quoted errors in figures represent only the statistical errors, which are obtained 
by a single elimination jack-knife method. Horizontal dashed lines represent the values in 
the infinite volume limit, which are evaluated in the case of n — 2, together with their one 
standard deviation. A summary of masses in the infinite volume limit, which are guided 
by various power-law behaviors, is given in Table IIVI In the case of the N* state, which 
receives the largest finite-size effect, the power two (n = 2) is most preferable because of 
X 2 /Ndf ~ 1 - 

For the nucleon case, we confirm that the finite-size effect is very small in the relatively 
heavier quark mass region. As shown in Fig. |T| all the data are located within lcr of the 
value in the infinite volume limit at k =0.1497 and 0.1506. We do not see any serious finite- 
size effect on the mass of the ground-state nucleon even at the smallest lattice size La ~ 1.6 
fm. In the mass of the negative-parity nucleon, a serious finite-size effect is not seen in the 
range of the spatial lattice size from 1.6 fm to 2.2 fm. This feature is consistent with that 
reported in Ref. jf|, in spite of the fact that their observed tendency of the size effect is 
opposite. Remarked that almost all subsequent lattice calculations to study the negative- 
parity nucleon are performed at the lattice size less than 2.2 fm. However, we find that the 
mass of the negative-parity nucleon suffers from the large finite-size effect in the range of 
the spatial lattice size from 2.2 fm to 3.2 fm. This means that the finite-size effect on the 
mass of the negative-parity baryon is unexpectedly severe and the spatial size is required 
at least 3 fm to remove this effect. What is the origin of the finite-size effect on the mass 
spectrum? In the phenomenological point of view, the “wave function” of baryon should be 
squeezed due to the small volume [l3(. The kinetic energy of internal quarks inside baryon is 
supposed to increase and thus the total energy of the three-quark system, which corresponds 
to the mass of baryon, should be pushed up. Such effect is expected to become serious for 
the excited state rather than the ground state. This intuitive picture seems to account for 
the observed trend of decreasing mass as the lattice size increases. 

In the spin-3/2 A channel, however, we find that the finite-size effect has a different 
pattern: the ground-state mass of the spin-3/2 A becomes large as the lattice size increases. 
This behavior may originate from a hyperhne interaction. This possibility will be discussed 
later. On the other hand, the negative-parity state of the spin-3/2 A, namely the Ag/ 2 , 
has the same pattern of the finite-size effect observed in the N* spectrum, while the finite- 
size correction of the Ag/ 2 is relatively smaller than that of the N*. As for the spin-1/2 
A channel, all data at different lattice sizes are roughly consistent with each other within 
relatively large errors. However, the significant finite-size effect may be hidden behind the 
large statistical errors. 

The finite-size effects on the masses of each baryon are summarized in Table 0 It is 
noteworthy that the finite-size corrections of A states ( J p = 3/2 ± ) and negative-parity 
nucleon N* at spatial size La= 1.6 fm can be seen even in the heavy quark region where 
that of the nucleon is almost negligible. Moreover, as we can learn from Figs RHTl the spatial 



lattice size is required to be as large as 3 fm to remove the finite-size effect. The systematic 
error stemming from the finite-size effect at spatial size La ~ 3.2 fm, ( (777.3.2 f m — 777 00 )/t 77 00 |, 
is smaller than 2% for all measured states. Therefore, in the later discussion, we analyze 
data obtained on the largest lattice size of 48 3 x 48 for spectroscopy of all hadrons. 

Finally, we would like to comment on other possible formula such as the exponential 
form: 

aMi = aMoo + cL” 1 exp(— L/L 0 ), (23) 


which is inspired by the Luscher’s formula for the asymptotic finite-size dependence of stable 


oy 

particle masses [12|. Here 1/Lo ~ M n can be expected in the phenomenological sense. The 
full three-parameter fits tend to yield considerably low values of 1 /L 0 . The obtained values 
are as low as ICC 4 in lattice units. It is clearly inconsistent with a relation 1 / L 0 rs./ M, 
since our simulations are performed in the range of M n ~ 0.2 — 0.5 in lattice units. If the 
parameter L 0 is fixed by 1/L 0 = M n , the resulting x 2 /^df from the two-parameter fits is no 
longer reasonable. Therefore, the finite-size behavior that we observed here can be described 
by a power-law formula El rather than above exponential formula. 


B. Chiral extrapolation 

All data of hadron masses computed on lattice with spatial size La ~ 3.2 fm are tabulated 
in Tables o and Em We perform a covariant single exponential fit to two-point functions 
of each hadron in respective fitting ranges. All fits have a confidence level larger than 0.15 
and x 2 /Ndf < 1-5. Our adopted ranges for fitting are basically determined by plateaus of 
each effective mass. 

Next, we extrapolate masses of all calculated hadrons to the chiral limit using two types 
of fitting formula: 


aM H — cq + c 2 (aM n )~ 

(24) 

( aMtf ) 2 = do T d 2 (nAAr) J , 

(25) 


where c 0 , c 2 , d 0 , d 2 are numerical constants. We evaluate the systematic error by the differ¬ 
ence between the chiral limit values obtained by two different fitting formulae. Hereafter, 
Eq. (12411 and (|25j) are referred as the “linear” fit and the “curve” fit, respectively. The chiral 
limit values and their values of fitting y 2 are listed in Table IVIIII For the p meson, the 
nucleon and the A 3 / 2 , in accordance with resulting y 2 values, we can determine that the 
curve fit is much preferable. However, for the other baryons, either fits yield acceptable 
values of y 2 because of the relatively large statistical errors on fitted data. However, the 
systematic error stemming from the difference of two fitting results is less than 10%. We 
then adopt the curve fit for all hadrons as the final analysis of chiral extrapolation. 

In Figs. EH we show squared masses of positive- and negative-parity baryons as a function 
of squared pion mass in each (/, J) channel. Circle symbols correspond to the negative- 
parity state (solid circles) and the positive-parity state (open circles). The extrapolated 
points in the chiral limit are represented by star symbols. In Table IIXI respective values in 
the chiral limit for all calculated baryons are listed and also compared with experimental 
values. Two types of input (r 0 input and M p input) are taken to expose the dependence 
on the choice of input to set a scale. We should keep in mind that the systematic errors 
stemming from this dependence are around 5%, the value of which exceeds the amount of 
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the statistical errors in the case of the nucleon and the A 3 / 2 . Instead, we quote various mass 
ratios, which do not suffer from such systematic uncertainties: 


Ma 3/2 /M n 

= 1.28(4)(9) 

(Expt. 

~ 1.31) 

Mn* / Mj\f 

= 1.61(10)(16) 

(Expt. 

~ 1.63) 

Ma* /2 / M n * 

= 1.28(7)(7) 

(Expt. 

~ 1.11) 


= 1.61(7)(13) 

(Expt. 

- 1.38) 


where the second quoted errors correspond to the systematic errors, which are estimated 
from the difference in the central values obtained by two types of chiral extrapolation. The 
mass ratio between the nucleon and its parity partner (or its hyperhne partner) shows a 
good agreement with the experimental values within statistical errors. However, mass ratios 
that include Ma* /2 are overestimated by about 15%. Our calculations are performed using 
the relatively heavy quark masses ( M n /M p ~ 0.66 — 0.96). The long chiral extrapolation is 
performed so that the evaluated values should not be taken too seriously. Indeed, both the 
“linear” fit and the “curve” fit do not include a term linear in aM n , which is responsible for 
the expected leading behavior close to the chiral limit in the quenched approximation. :27] 
Finally we stress that the level ordering in A spectra 1 , Ma 3/2 < Ma* /2 ~ Ma* /2 < Ma 1/2 , 
is well reproduced in comparison to corresponding A states, which are all ranked as four 
stars on the Particle Data Table j2,gj. In addition, it is worth mentioning that a signal of 
A(1750) (/ = 3/2 and J p = l/2 + ), which is the weakly established state (one star) |28[, 
cannot be seen in our data. 


C. Hyperfine mass splitting 

In this subsection, we discuss our numerical results on the hyperfine mass splittings 
in both parity channels ( e.g. Ma 3/2 — M N and Ma* /2 — M N *). In the quark potential 
model, if the inter-quarks potential is central force and independent of the flavor and spin, 
the hamiltonian has spin-flavor SU(6) symmetry |29[. Under the exact spin-flavor SU(6 ) 
symmetry, masses of resonance in the same SU(6) multiplet, e.g. 56-plet for the nucleon and 
the A 3 / 2 particle or 70-plet for negative-parity N* and A* , 2 states, should be degenerate. 
However, the actual baryon spectra show evident violations of this symmetry. In the non- 
relativistic quark models, violations of the spin-flavor SU(6) symmetry are caused by spin- 
dependent interactions [291 ]. 

First, we consider the finite-size effect on the hyperhne mass splittings (N — A 3 / 2 and 
N* — Ay 2 ). The lattice-size dependences of the hyperhne mass splittings are shown in Figs. [7] 
and E We observe an unique feature in both parity channels. Each hyperhne splitting 
becomes small as the spatial lattice size decreases 1 2 . This feature is rather peculiar from the 
viewpoint of the non-relativistic quark models. The hyperhne interaction may be derived 
from one-gluon exchange as a spin-spin component of the Fermi-Breit type interaction. 
Thus, the hyperhne mass splitting originates from a contact interaction between colored 


1 In our results, AU 2 and A*/ 2 , both of which belong to the same SU( 6) multiplet, are mostly degenerate 
within statistical errors. Strictly speaking, we have not seen clear splitting between them. 

2 The same feature is observed in the charmonium spectrum l.'Ifl ] . The hyperfine splitting between g c and 
J/ip diminishes as the spatial lattice size decreases if La Ss 1.3 fm. 
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quarks @]. As the size of a hadron decreases, probability of finding two quarks at the same 
spatial point inside hadron increases, so that the hyperfine mass splitting would increase. 
However, the observed finite-size effect where the hyperfine mass splitting diminishes as the 
lattice size decreases is opposite from above naive expectation. At present, we do not have a 
simple interpretation of the finite-size behavior of the hyperfine mass splitting in our data of 
quenched lattice QCD. This peculiar behavior may suggest some other origin of the hyperfine 
interaction rather than one-gluon exchange. Finally, it is worth mentioning that the 
observed size dependence of the hyperfine mass splitting accounts for the opposite pattern 
of the finite-size effect on masses of the A 3 / 2 and also our observation that the finite-size 
effect for the A* , 2 state is slightly milder than the N* state. 

We shortly comment on the splittings between pairs of parity partner (N — N* and 
A3/2 — Ag, 2 ). It is found that these splittings become small as the lattice size decreases 
in both N and A channels. All data including hyperfine mass splittings are tabulated in 
Table H 

Next, we discuss the quark-mass dependence of the hyperfine mass splitting. In Fig. 0 
the hyperfine mass splittings between the nucleon and the A 3 / 2 are plotted as a function of 
sum of the nucleon and A 3 / 2 masses by using data at the largest spatial lattice size, La ~ 3.2 
fm. We find that the hyperfine mass splitting becomes large as the quark mass decreases. 
This feature shows a good agreement with the general form of the spin-dependent interaction 
in the non-relativistic quark models where the spin-dependent interaction is proportional to 
inverse powers of the quark mass 2Q\. Of course, one can easily find such the feature in the 
actual baryon spectra in the wide range from the light (up, down) sector to the charm sector. 
Fig. El includes some experimental points (stars), which correspond to spin 1/2 and 3/2 
doublets; A(939)-A(1232), £(1192)-£(1385), £(1315)-£(1530) and £ c (2455)-£ c (2520). 
Our data including the chiral extrapolated value, Ma 3/2 — M N = 0.262(11) GeV, is fairly 
consistent with those experimental points. It is also found that the quark-mass dependence 
for the hyperfine mass splitting between the Ag , 2 and the N* is consistent with the naive 
expectation, which is deduced from the explicit mass dependence of the spin-dependent 
interaction in the non-relativistic quark models. 

Finally, it is important to note that the value of the hyperfine mass splitting should be 
sensible to the leading discretization errors of the Wilson fermion action, which may induce 
the extra chromomagnetic-moment interaction at finite lattice spacing. We however believe 
that qualitative features of the finite-size dependence and the quark-mass dependence, which 
are observed in this paper, do not change with respect to discretization 0(a) errors. 


IV. CONCLUSIONS 

In this paper, we have studied the finite-size effect on masses of various baryons: the 
nucleon states in both parity channels ( J p = l/2 ± ) and the A states in all spin-parity 
channels ( J p = l/2 ± ,3/2 ± ). Our quenched lattice simulations were employed at relatively 
weaker coupling, (3 = Q/g 2 = 6.2, where the cut-off scale (~ 3 GeV) is definitely higher than 
mass scale of all observed baryons (~ 1 — 2 GeV). Three different lattice sizes, La ~ 1.6, 2.2 
and 3.2 fm, were utilized to examine how large lattice size is required for excited baryon 
spectroscopy. 

We have found the considerable finite-size effect on masses of all A states and negative- 
parity nucleon N* even in the relatively heavy quark region where the finite-size effect on the 
nucleon is almost negligible. The finite-size behavior that we observed can be described by a 
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power law M^ — Ml oc L~ n with n ~ 2 rather than the ex pon ential form oc L exp(— L/L 0 ). 
This observation is consistent with that reported in Ref jl3(. If the f ini te-size effect is kept 
as small as a few percent level, the spatial lattice size La ^ 3 fin should be required for 
excited baryons, especially for the negative-parity nucleon. 

The finite-size behavior of the power law might originate from the phenomenological point 
of view as the squeezed hadron size. Indeed, we confirmed that the rather large finite-size 
effect appears for excited baryons, the “wave function” of which ought to be extended rather 
than that of the ground state. The squeezed “wave function” may increase the total energy 
of the three-quark system because of a gain in kinetic energy of internal quarks inside hadron 
so that the mass is expected to be pushed up as the lattice size decreases. However, the 
A3/2 state, which is the lowest lying state in the A channel, reveals the opposite pattern of 
the size effect. The mass of the A 3 / 2 decreases as the lattice size decreases. This peculiar 
behavior might attribute to the finite-size effect on the hyperhne interaction. 

According to the SU(6) quark model, the nucleon and the A 3 / 2 belong to the same SU(6 ) 
multiplet so that the difference of wave functions between the N and the A 3 / 2 is induced by 
the hyperhne interaction. Our study of the finite-size effect were performed in the relatively 
heavy quark mass region, where the finite-size effect on the nucleon mass is almost negligible. 
Therefore, the finite-size effects on the A 3 / 2 can attribute to the effect of the hyperhne mass 
splitting. Indeed, in the case of the N* — A^ 2 hyperhne mass splitting, we observed the same 
pattern of the finite-size effect where N* — A* , 2 splitting diminishes as the lattice spatial 
size decreases, while either N* or A^ states yield the expected size effect where both 
masses increase as the lattice size decreases. However, the hyperhne interaction induced 
by one-gluon exchange, which is mainly a contact interaction between colored quarks, may 
not correctly account for this peculiar behavior of the finite-size effect. It seems that our 
observed finite-size effect on the hyperhne mass splitting suggests some other origin of the 
hyperhne interaction for baryons rather than one-gluon exchange. 
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APPENDIX A: (ANTI-)PERIODIC TIME DEPENDENCE FOR FINITE T 

Let us consider the functional form of the typical two-point correlation function G(t) = 
e —M|t| unf ] er the periodic (anti-periodic) boundary condition in time. The desired function, 
in the finite extent T with the periodic boundary condition, should be formed as 

OO 

Gp.h.c.if) = G(t + nT), 


n =—oo 


(Al) 
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which certainly satisfies G p .b.c.(f) = G p ^_ c _(t + T). Recall that the given correlation function 
G(t) clearly has the time-reflection symmetry G(t) = G(-t) 1 . Thus, one can rewrite G p ^ c ft) 
for 0 < t < T as 


<W<) = V G(t + nT) + G ( nT - t) 


n =0 
oo 

E' 

71=0 

OO 

E 

n =0 


-M{t+nT) 


n=l 

oo 

E' 

n= 1 


-M(nT-t) 


e -nMT j-g-Mt + e -M(T-t)j 


(A2) 


In the case of the anti-periodic boundary condition, where G a . P .b.c.(£) = — G a .p.b.c.(^ + T), 
one obtains the following general form: 


Ga.p.b.e.00 = X] (~TG(t + 7lT) 

n =— oo 

oo oo 

= £(—)”G(t + nT) + J2(-) n G(nT - t) 


71=0 

OO 


71=1 


n —tiMT f —Mt _ —M(T—t 




71=0 


(A3) 


We would like to make the following few comments. First, we can derive familiar forms 
from Eqs. M) and m as the asymptotic form in the large T limit (MT >> 1). 


G 


p.b.c./a.p.b.c 


. 0 ) 


-Mt 


± g-M(T-t) 


which has only a contribution from the first wrap-round effect. The (n — l)-th wrap-round 
effect should be suppressed with the factor e~ nMT . Then, the linear combination of G p ±,. c ft) 
and G a . p .b.c.(f) yields 

G(t) = — [G p .b.c. {t) + G a . p .b.c.(f)] ~ e Mt 

where the first wrap-round effect is canceled. 

Secondary, we will see that G(t) corresponds to a periodic function in the finite extent 
2 T as follows. One sums up all wrap-round contributions in Eqs. m and <ra, and then 
obtains the exact functional forms: 


Gp.b.cXO 

G a.p.b.c. if) 


1 r -Mt , —M(T—t)~\ _ cosh [M(t 2 )] 

l-e~ MT[ J sinh[M|] 

1 r -Mt _ _ sinh[M(f — -y)] 

1 + e~ MT V - J cosh[M|] 


(A4) 

(AS) 


1 In the case of G(t) = sgn (t)e M \ t \ J a property G(t) = — G(—t ) should be utilized instead of the time- 
reflection symmetry. 
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where the (anti-)periodic time dependence for finite T is given by a cosh(sinh) function. 
Therefore, G(t) is given by 

77/^ _ 1 ( cosh [M(t - f)] sinh[M(f - f )] \ _ cosh[M(f - T)] 

2 y sinh[M-|] cosh[M^] J sinh[MT] 

The hnal expression is nothing but G p , b. c .(t) with the 2T-periodicity. 
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TABLE I: Simulation parameters for each volume studied in this work. 


p 

lattice size (L 3 x T ) 

kappa values 

statistics 

~ La (fm) 

6.2 

24 3 x 48 

{0.1506, 0.1497} 

320 

1.6 


32 s x 48 

{0.1506, 0.1497} 

240 

2.2 


48 3 x 48 

{0.1520, 0.1506, 0.1497, 0.1489, 0.1480} 

210 

3.2 


TABLE II: The physical scale set by the Sommer parameter ro = 0.5 fm. The value of r^/a is 
taken from Ref. Eh], 


— 

P 

r 0 /a 

oT 1 (GeV) 

a (fm) 

6.2 

7.380(3) 

2.913 

0.06775 

.. ..... ...... ..... .. —..... ' - -- 



TABLE III: Finite-size 

effect on masses of all measured baryons at /3 = 6.2 


L 3 x T 

K 

aAfjv 

a Mn* 

aM A 3/2 

aM A * /2 

aM A 1/2 

aM A * /2 

24 3 x 48 

0.1506 

0.1497 

0.594(5) 

0.667(4) 

0.829(19) 

0.891(15) 

0.624(7) 

0.693(5) 

0.872(16) 

0.928(13) 

0.993(24) 

1.055(22) 

0.777(28) 

0.864(22) 

32 s x 48 

0.1506 

0.1497 

0.597(4) 

0.670(3) 

0.808(15) 

0.873(13) 

0.628(7) 

0.697(4) 

0.878(13) 

0.928(12) 

0.985(27) 

1.043(25) 

0.835(25) 

0.901(21) 

48 3 x 48 

0.1506 

0.1497 

0.594(4) 

0.668(3) 

0.762(15) 

0.832(13) 

0.637(5) 

0.702(4) 

0.837(14) 

0.892(12) 

0.964(26) 

1.023(25) 

0.814(28) 

0.869(23) 
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TABLE IV: Infinite volume limit of various baryon masses at (3 = 6.2 


K 

type 

aM N 

X 2 /N D f 

aM]y* 

X 2 /Ndf 

aM A y 2 

X 2 /Ndf 

aM/\* 

3/2 

X 2 /Ndf 

0.1506 

L~ 3 

0.597(4) 

0.46 

0.763(16) 

1.54 

0.638(5) 

0.42 

0.846(14) 

3.06 


L~ 2 

0.595(5) 

0.47 

0.747(20) 

0.94 

0.641(7) 

0.26 

0.836(18) 

2.56 


L~ l 

0.595(9) 

0.47 

0.697(36) 

0.45 

0.651(12) 

0.12 

0.808(31) 

2.01 

0.1497 

L~ 3 

0.669(4) 

mm 

0.833(13) 

1.70 

0.703(4) 

0.22 

0.898(12) 

2.51 


L~ 2 

0.669(5) 

mt 

0.819(17) 

1.06 

0.705(6) 

0.11 

0.889(16) 

2.00 


L- 1 

0.669(8) 

BSS 

0.775(30) 

0.52 

0.712(10) 

0.03 

0.862(27) 

1.47 


TABLE V: Summary table of the finite-size effect on various baryon masses for k = 0.1506. 


baryon 

as increase L 

~ (Afi.efm — Af 00 )/M 00 

~ IM3.2 frn — Moo |/Moo 

N 

-> 

less than 1% 

less than 1% 

N* 

\ 

11% 

2% 

A3/2 

/ 

-3% 

less than 1% 

A* 

^ 3/2 

\ 

4% 

less than 1% 


TABLE VII: Fitted 

masses of A baryons ( J p 

= 3/2 ± 

, l/2 ± ) at 0 = 

6.2 and lattice volume 48 3 x 48. 

K 

o.Ma 3 / 2 

range 



«Ma 1/2 

range 

CL Ad A* 

^1/2 

range 

K c 

0.424(11) 

N/A 

0.681(25) 

N/A 

0.798(38) 

N/A 

0.640(53) 

N/A 

0.1520 

0.535(9) 

[22,28] 

0.770(21) 

[12,20] 

0.883(34) 

[10,15] 

0.732(54) 

[13,20] 

0.1506 

0.637(5) 

[22,28] 

0.837(14) 

[12,20] 

0.964(26) 

[10,15] 

0.814(28) 

[13,20] 

0.1497 

0.702(4) 

[22,28] 

0.892(12) 

[12,20] 

1.023(25) 

[10,15] 

0.869(23) 

[13,20] 

0.1489 

0.759(4) 

[22,28] 

0.943(12) 

[12,20] 

1.076(25) 

[10,15] 

0.921(20) 

[13,20] 

0.1480 

0.822(4) 

[22,28] 

1.001(11) 

[12,20] 

1.136(25) 

[10,15] 

0.981(18) 

[13,20] 


TABLE VI: Fitted masses of pion, p -meson and nucleon ( J p = l/2 ± ) states at fi = 6.2 and lattice 
volume 48 3 x 48. 


K 

aMjr 

range 

aM p 

range 

aM]\f 

range 

aMjy* 

range 

K c 

0 

N/A 

0.248(3) 

N/A 

0.332(9) 

N/A 

0.534(33) 

N/A 

0.1520 

0.2141(8) 

[15,30] 

0.324(2) 

[14,30] 

0.469(5) 

[20,28] 

0.640(22) 

[13,20] 

0.1506 

0.3167(9) 

[15,30] 

0.388(1) 

[14,30] 

0.594(4) 

[20,28] 

0.762(15) 

[13,20] 

0.1497 

0.3726(9) 

[15,30] 

0.430(1) 

[14,30] 

0.668(3) 

[20,28] 

0.832(13) 

[13,20] 

0.1489 

0.4188(9) 

[15,30] 

0.467(1) 

[14,30] 

0.730(3) 

[20,28] 

0.891(12) 

[13,20] 

0.1480 

0.4678(9) 

[15,30] 

0.509(1) 

[14,30] 

0.799(3) 

[20,28] 

0.957(11) 

[13,20] 
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TABLE VIII: Results of the chiral extrapolated masses in lattice units for all measured hadrons 
by two types of fitting formula, the linear and curve fits. 


type 

aM p 

aMjv 

uMn* 

aM A3/2 

aM A * , 

3/2 

aM Al/2 

aM /\* 

^ 1/2 

linear 

(mass) 

0.283(2) 

0.403(7) 

0.582(26) 

0.477(8) 

0.703(22) 

0.818(35) 

0.670(44) 


( x 2 /Ndf) 

6.20(1.67) 

9.62(3.37) 

0.62(0.45) 

1.57(0.99) 

0.03(0.08) 

0.01(0.02) 

0.01(0.05) 

curve 

(mass) 

0.248(3) 

0.332(9) 

0.534(33) 

0.424(11) 

0.681(25) 

0.798(38) 

0.640(53) 


(. X 2 /Ndf ) 

1.59(0.63) 

0.03(0.14) 

0.07(0.12) 

0.05(0.15) 

0.17(0.21) 

0.01(0.03) 

0.01(0.05) 


TABLE IX: The third and forth columns list results of all measured baryon masses in GeV units, 
which are set by two different inputs, ro input and M p input. The fifth and sixth list experimental 
values of the corresponding baryon and its status in the Particle Data Table [3]. The possible 
assignments of the SU{ 6) <g) 0(3) supermultiplet are also embedded into the final column. 


state 

(I: J P ) 

our results [GeV] 


physical baryon 

(status) 

517(6) ® 0(3) 



(r 0 input) 

(M p input) 



classification 

N 

~7TT T 7 

V2> 2 > 

0.967(26) 

1.030(27) 

N( 939) 


[56,0+] 

N* 

(1 

V 2 ’ 2 > 

1.555(96) 

1.658(102) 

JV(1535) S n 


[70,1-] 

A3/2 

(A 3 + ) 
i 2 ’ 2 > 

1.235(32) 

1.316(34) 

A(1232) P 33 


[56,0+] 

A* 

^ 3/2 

(3 3 —\ 
i 2 ’ 2 > 

1.983(72) 

2.114(77) 

A(1700) D 33 


[70,1-] 

Al/2 

(3 1+) 

2.324(110) 

2.477(117) 

A(1910) P 31 


[56,2+] or [70,0+] 

A* 

^1/2 

(3 I") 

' 2 ’ 2 ) 

1.864(154) 

1.987(164) 

A(1620) S 31 


[70,1-] 


TABLE X: Finite-size effect on various mass splittings at /3 = 6.2 


L 3 xT 

Av 

aMA 3/2 -aMjv 

aM A * /2 -aMN* 


aM A * /2 -aM A3/2 

aM Al/2 -aM A * /2 

24 3 x 48 

0.1506 

0.1497 

0.030(7) 

0.026(5) 

0.043(20) 

0.037(14) 

0.235(21) 

0.224(16) 

0.248(18) 

0.235(15) 

0.216(36) 

0.191(30) 

32 s x 48 

0.1506 

0.1497 

0.030(6) 

0.027(4) 

0.070(17) 

0.055(13) 

0.211(17) 

0.201(14) 

0.250(15) 

0.231(13) 

0.150(37) 

0.143(33) 

48 3 x 48 

0.1506 

0.1497 

0.043(5) 

0.035(4) 

0.075(16) 

0.060(12) 

0.168(18) 

0.164(15) 

0.200(16) 

0.190(14) 

0.151(35) 

0.154(31) 
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FIG. 1: Masses of nucleons ( J p = l/2 ± ) in lattice unit as functions of spatial lattice size in the 
physical unit for k = 0.1506 (left figure) and k = 0.1497 (right figure). Solid curves are fits of the 
form clMl, = aM^ + cL~ n with the value n = 2. Horizontal dashed lines represent extrapolated 
values for each parity state in the infinite volume limit, together with their one standard deviation. 




FIG. 2: Masses of spin-3/2 A baryons in lattice unit as functions of spatial lattice size in the 
physical unit for k = 0.1506 (left figure) and n = 0.1497 (right figure). Solid curves and dashed 
lines are defined as in Figured 
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FIG. 3: Masses of spin-1/2 A baryons in lattice unit as functions of spatial lattice size in the 
physical unit for k = 0.1506 (left figure) and k = 0.1497 (right figure). Solid curves and dashed 
lines are defined as in Figured 



(aMJ 2 


FIG. 4: Squared masses of nucleons ( J p = 1/2 1 * 1 ) in lattice unit as functions of squared pion mass. 
Circle symbols correspond to the negative-parity state (solid circles) and the positive-parity state 
(open circles). The extrapolated points in the chiral limit are represented by stars respectively. 
The experimental values are marked with lower and upper open diamonds. 
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FIG. 5: Squared masses of A baryons ( J p = 3/2 ± ) in lattice unit as 
mass. All symbols are defined as in Figure El 



FIG. 6: Squared masses of A baryons ( J p = l/2 ± ) in lattice unit as 
mass. All symbols are defined as in Figure El 


functions of squared pion 


functions of squared pion 
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FIG. 7: Masses of spin-1/2 A baryons in lattice unit as functions of spatial lattice size in the 
physical unit for k = 0.1506 (left figure) and k = 0.1497 (right figure). Solid curves are fits of the 
form ciMl = aM, x + cL~ n with the value n = 2. Horizontal dashed lines represent the mean value 
of the infinite volume extrapolation and one standard deviation from it. 




FIG. 8: Masses of spin-1/2 A baryons in lattice unit as functions of spatial lattice size in the 
physical unit for k = 0.1506 (left figure) and k = 0.1497 (right figure). Solid curves and dashed 
lines are defined as in Figured 
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aM 4 +aM N 


FIG. 9: Hyperfine splitting between the N and A 3 / 2 states as a function of Ma + Mjy. All 
calculations are done for three degenerate valence-quarks. The experimental values in the wide 
range from the light (up, down) sector to the charm sector are marked by stars. 




